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Abstract 

Numerical  results  demonstrate  that  both  Fraunhofer  and  Fresnel  diffractions  provide 
good  approximations  to  forward  scattering  by  elastic  spherical  shells  and  rigid/soft 
spheres,  for  non-dimensional  frequency  values  5  <  ka  <  40,  and  scattering  angles 
0°<6<  J0°.  The  calculations  were  done  for  the  far  field,  with  incident  plane  waves. 


Introduction 


Forward  scattering  has  been  used  in  particle  sizing  in  optical  [I]  and  acoustic  [2] 
measurements.  Both  Fraunhofer  diffraction  and  Fresnel  diffraction  [3]  have  been  used  as 
predictive  tools  in  calculations  of  scattering  from  objects  of  suitably  small  sizes.  Since 
acoustic  waves  penetrate  elastic  objects,  but  the  Fraunhofer  and  Fresnel  diffractions  do  not 
include  elasticity,  the  validity  of  these  physical  optics  approaches,  even  in  forward 
scattering,  may  be  questioned.  In  this  report  we  compare  predictions  of  forward  scattering 
from  elastic  spherical  shells  and  rigid/soft  spheres,  using  known  exact  formulas  [4],  to 
Fraunhofer  and  Fresnel  predictions. 

The  methods  of  physical  optics  have  been  'iplied  to  acoustic  scattering  problems  by 
many  authors  [2,  5-12]  with  varying  degrees  of  success.  The  best  known  of  these 
techniques  is  the  Kirchhoff  approximation  [5]  which  assumes  that  the  incident  field  and  its 
normal  derivative  are  unperturbed  in  an  open  aperture,  and  are  zero  on  the  shadow  side  of 
a  scatterer.  The  application  of  these  assumptions  to  backscattering  from  semi-infinite  plates 
[5],  and  from  spheres  and  cylinders  [6-8]  has  shown  that  the  Kirchhoff  method  is 
unreliable  in  these  instances. 

Application  of  the  Kirchhoff  method  to  forward  scattering  from  a  semi-infinite  plate  has 
been  discussed  by  Officer  [9].  Pierce  and  Hadden  find  [10,  11]  that  for  small  forward 
angles,  the  Kirchhoff  method  provides  a  good  approximation  to  scattering  by  a  wedge. 
Hunter,  Lee,  and  Waag  [2]  have  experimentally  measured  the  forward  scattering  patterns  of 
single  nylon  filaments.  They  demonstrated  that  the  patterns  were  well  explained  by  a 
Huygens  construction,  calculated  numerically.  The  Fraunhofer  and  Fresnel  calculations  [3] 
discussed  below  are  subsets  of  the  Huygens  construction.  Kirchhoff  predictions  for  forward 
scattering  by  an  infinite  rigid  cylinder  have  been  shown  by  one  of  the  authors  [12]  to  agree 
with  exact  predictions  [13].  For  forward  scattering,  the  Kirchhoff  method  is  identical  to 
Fresnel  diffraction  [3,  9].  Since  both  Fresnel  and  Fraunhofer  diffractions  are  special  cases 
of  a  unified  theory  based  on  the  Fresnel-Kirchhoff  integral  [3],  we  prefer  to  retain  the 
Fresnel  label  instead  of  the  Kirchhoff  label.  We  show  that  in  the  examples  considered  here, 
both  Fraunhofer  and  Fresnel  diffractions  provide  good  approximations  to  forward 
scattering  by  elastic  spherical  shells  and  rigid/sofl  spheres  [14]. 

1.  Exact  formulas 

The  well-known  exact  formulas  [4]  for  spherical  elastic  and  rigid/soft  scattcrers  arc 
summarized  in  the  present  section.  For  spherical  scattcrers,  the  complex  pressure  of  the 
incident  plane  wave  can  be  written  as 

Pi^.  =  expfikrcosB),  (1) 


■ion  for 

CRA&I 

TAB 

’  C.-l’.iC'T' 


where  k  is  the  wave  number,  r  is  the  distance  of  the  observation  point  from  the  center  of  ^  ^  .  |  .hii  ly  '^'odes 
the  sphere,  and  0  is  the  spherical  polar  angle,  '.vith  0  =  0  for  the  forward  direction.  ^  _ _ 
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Similarly,  for  the  scattered  signal 

Pscm  =  “  (^)®*P  +  l)A,P,(cosfl)  =  (fjexp  (ilcr)f(Q)  .  (2) 

Here  P„(cos6)  are  Legendre  polynomials,  and  the  coefficients  A„  are  determined  from  the 
elastic  boundary  conditions  [4|.  The  formulas  are  for  the  asymptotic  case  kr  >  I.  For  rigid/ 
soft  boundary  conditions  [4] 

A„(rigid)  =  -j;(ka)/h'„(ka),  and  A„(soft)  =  -j„(ka)/h„(ka) .  (3) 


Here  the  are  spherical  Bessel  functions,  and  the  h„  are  spherical  Hankel  functions  of  the 
first  kind.  The  total  field  at  any  point  is  given  1^  the  coherent  sum  of  the  incident  and 


scattered  fields  as 


Pu*  Pw  Pscao  • 


(4) 


2.  Physical  optics  formulas 

In  the  physical  opt'cs  approach  to  the  forward  scattering  problem,  we  replace  the  scattering 
object  (sphere)  with  an  opaque  screen  having  the  same  area  of  cross  section.  Consider  the 
case  when  the  wavefront  due  to  the  incident  signal  coincides  with  the  screen.  The  complex 
pressure  at  any  forward  point  is  then  given  by  the  Fresnel-Kirchhoff  integral  [3] 


Pw  = //expIikf(^,Tl)d^dTil 


(5) 


for  incident  plane  waves  of  unit  strength  and  small  forward  scattering  angles.  Here  k  is  the 
wave  number,  and  denote  the  two-dimensional  cartesian  co-ordinates  of  a  point  Q  on 
the  opaque  screen  (Figure  I).  Let  the  point  P'  denote  the  source  (at  infinity),  the  point  P 
denote  the  receiver  (observation  point),  and  the  point  O  denote  the  intersection  of  the 
straight  line  P'P  with  the  opaque  screen.  Let  s'  =  OP,  and  s  =  QP.  Then  the  function 
f(4.h)  =  s  -  s'.  The  integral  in  Eq.(5)  runs  over  the  entire  insonified  region,  viz.  the 
wavefront  extending  to  too  in  and  q,  but  excluding  the  opaque  screen. 


Figure  1. 

Source,  receiver,  co-ordinate 
axes,  and  the  position  of  the 
wavefront  as  it  coincides 
with  the  screen 


For  better  physical  insight,  let  us  rewrite  Eq.(S)  as 


Pw  =  C'J  f exp  likf(^,q)ld^dq  -  C'jjcxp  [ikf(^,q)Jd^dq  =  -  pji^  ,  (6) 

where  the  integral  labeled  ‘disc’  runs  over  the  opaque  screen.  Equation  (6)  can  be 
understood  as  a  result  of  Babinefs  principle  [3].  The  first  term  pi„^.  in  Eq.(6)  represents  the 
field  due  to  the  entire  infinite  wavefront,  i.e.,  the  field  in  the  absence  of  any  scatterer.  The 
second  term  pjj^  in  Eq.(6)  represents  the  field  due  to  the  complementary  geometry,  i.e.,  an 
insonified  opening  in  an  infinite  opaque  screen.  The  terms  p,,^  and  pj„^.  in  Eq.(6)  have 
exactly  the  same  meaning  as  the  corresponding  quantities  in  Eq.(4). 
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Within  the  Fraunhofer  upproximations.  using  some  change  of  variables  and 
straightforward  int^ration,  we  obtain  13] 

_  /c  I.  r  .  i  fexpOks')!,,  ,,r2J|(kasine)] 
p,.,(Fraunhofer)  =  -  -  [— j^J  ]  •  (7) 

Here  6  is  the  angular  separation  of  a  point  on  the  diffraction  pattern  measured  from  the 
forward  direction,  and  is  the  same  as  the  angle  6  in  the  previous  section. 

To  calculate  p^j^^  using  the  Fresnel  approximation.  Bom  and  Wolf  [3]  choose  a 
co-ordinate  system  such  that  the  linear  terns  in  f(^,Ti)  are  identically  zero.  Defining  two 
new  variables  u  =  ^(k/ns')''^,  and  v  =  T\(k/Rs')'^,  we  can  write  [3] 

PdijcCFresnel)  =  -  y  exp  (iks')tCu,(disc)  +  iS„^(disc)l  ,  (8) 

where 

C„„(disc)  =  2 J  cos  C lv^(u)]du  -  2  Jsin  [-j  u^)  S[v^(u)]du  (9) 

due  disc 

S„v(disc)  =  if  sin  (-1  u^)  C  Iv^(u)]du  +  if  cos  (-1  u^)  S[v^(u)]du  (10) 

diac  disc 

v™a*  =  [a'^  -  (u  +  Uof  1“  .  a'  =  [-Kt:  +  7)]  a  (11) 

and  Uq  is  the  (positive)  distance  in  the  u-v  plane  between  the  origin  of  the  co-ordinate 
system  and  the  center  of  the  opaque  disc.  This  distance  is  also  equal  to  the  distance  in  the 
observation  plane  (also  measured  in  u-v  space)  of  a  point  on  the  diffraction  pattern  from 
the  center  of  the  pattern.  The  quantities  C  and  S  on  the  right  hand  sides  of  Eqs.(9)  and  (10) 
are  the  well  known  Fresnel  integrals  [3,  IS].  We  calculated  C„v(disc)  and  S„.(disc)  of 
Eqs.(8)-(10)  by  numerical  integration  using  Simpson's  rule. 

The  incident  field  pj^^  is  a  plane  wave  of  unit  strength,  given  l^ 

p,^  (Fraunhofer)  =  exp  (iks'eosO)  (*2) 

in  polar  co-ordinates  (s',6),  used  in  the  Fraunhofer  calculation.  As  mentioned  above,  the 
Fresnel  calculation  uses  a  special  co-ordinate  system  [3]  in  which  the  diffracted  field  is 
calculated  on  a  plane  perpendicular  to  the  incident  signal  direction,  at  a  distance  s'  from 
the  scatterer.  For  this,  the  incident  signal  is  given  1^ 

Pi,^, (Fresnel)  =  exp  (iks')  .  (13) 


3.  Numerical  results 

Here  we  present  the  results  of  our  calculation  of  the  total  field  based  on  the  Fraunhofer  and 
the  Fresnel  approximations.  These  results  are  compared  with  corresponding  results 
obtained  using  exact  formulas  [4]  for  spherical  steel  shells  and  rigid/soft  spheres.  In  all 
examples,  we  have  calculated  the  diffraaion  pattern  in  the  fer  field,  with  (s'/X)  =  I(XX),  for 
incident  plane  waves.  In  the  Fresnel  calculation,  the  increment  for  numerical  integration  in 
u-space  was  set  at  Au  =  0.001;  the  fields  were  calculated  at  u-space  intervals  of  Au  =  0.04. 
This  generated  198  plotting  points  in  the  angular  range  0°  <  0  <  10°.  In  the  Fraunhofer  and 
exact  elastic/rigid/sofi  calculations,  the  fields  were  calculated  at  angular  intervals  of 
A0  =  0.05°.  For  the  external  fiuid  medium,  water,  sound  speed  =  1482.5m/s;  density  = 
IO(X)kg/m^.  For  steel,  longitudinal  sound  speed  =  5950m/s;  shear  speed  =  3240m/s; 
density  =  77(X)kg/m^.  For  air  inside  the  steel  shell,  sound  speed  =  330m/s;  density  = 
1.3kg/m3. 

Figures  2—5  show  the  squared  magnitude  of  the  total  field  p,^  plotted  vs.  the  forward 
scattering  angle  0,  for  ka  =  5  and  40.  Since  signal  penetration  into  the  scatterer  would  be 
greater  at  larger  wavelengths,  the  choice  of  ka  =  5  (wavelength  =  1.26  x  sphere  radius)  is 
intended  to  demonstrate  that  the  physical  optics  predictions  are  valid  even  at  such  relatively 
large  wavelengths.  The  Fraunhofer  and  Fresnel  predictions  arc  compared  with  rigid/soft 
results  in  Figures  2  and  3.  Exact  predictions  for  steel  shells  of  1%,  5%,  and  10% 
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uhicknesses  are  compared  with  the  Fresnel  results  in  Figs.  4  and  S.  Angular  positions  of  the 
diffraction  crests  and  troughs  predicted  by  the  four  approaches,  viz.  Fraunhofer,  Fresnel, 
steel  shell,  and  rigid/soft  show  good  mutual  agreement  in  the  angular  positions  of  the  crests 
and  troughs.  In  all  cases  the  a^mptotic  value  (as  6  increases)  is  OdB,  consistent  with  the 
normalization  of  the  incident  signal  strength  to  unity,  and  the  fact  that  the  scattered  field 
strength  &lls  off  rapidly  for  increasing  6.  The  unit  normalization  of  the  incident  signal  also 
explains  why  the  oscillations  are  centered  at  OdB. 


Fraunhofer 


Figure  2.  >Mation  of  the  total  acoustic  field  strength  with  forward  scattering  angle  6  for  rigid/solii 
spheres,  compared  with  Fraunhofer  and  Fresnel  diffractions,  ka  =  S. 


I  1 


i  -0.5- 


-  Fraunhofer 

-  Fresnel 

.  rigid 
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0  (deg) 


Figure  3.  Similar  to  Figure  2,  for  ka  -  40. 
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Figure  4.  \^iiatian  of  the  total  acoustic  field  strength  with  forward  scattering  angle  6  for  1%,  5%  ahd 
10%  steel  shells  compared  with  Fresnel  diffraction,  ka  =  S. 


Figure  5.  Similar  to  Figure  4,  for  ka  =  40. 


4.  Conclusions 

Both  Fraunhofer  and  Fresnel  diffractions  have  been  demonstrated  to  provide  good 
approximations  to  forward  scattering  by  elastic  spherical  shells  and  rigid/soft  spheres,  for 
non-dimensional  frequency  values  5  <  ka  <  40,  and  scattering  angles  0°  <  8  <  10°. 
Calculations  were  done  for  the  far  Held,  for  incident  plane  waves.  Angular  positions  of  the 
diffraction  crests  and  troughs  predicted  by  the  four  approaches,  viz.  Fraunhofer.  Fresnel, 
elastic  shell,  and  rigid/soft  arc  in  good  mutual  agreement. 
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